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The Extragalactic Background Light (EBL) is the integrated light from all the stars that have 
ever formed, and spans the IR-UV range. The interaction of very-high-energy (VHE: E > 100 GeV) 
7-rays, emitted by sources located at cosmological distances, with the intervening EBL results in 
e~e + pair production that leads to energy-dependent attenuation of the observed VHE flux. This 
introduces a fundamental ambiguity in the interpretation of the measured VHE blazar spectra: 
neither the intrinsic spectra, nor the EBL, are separately known - only their combination is. In 
this paper we propose a method to measure the EBL photon number density. It relies on using 
simultaneous observations of blazars in the optical, X-ray, high-energy (HE: E > 100 MeV) 7-ray 
(from the Fermi telescope), and VHE 7-ray (from Cherenkov telescopes) bands. For each source, 
the method involves best-fitting the spectral energy distribution (SED) from optical through HE 7- 
rays (the latter being largely unaffected by EBL attenuation as long as 2^ 1) with a Synchrotron 
Self-Compton (SSC) model. We extrapolate such best-fitting models into the VHE regime, and 
assume they represent the blazars' intrinsic emission. Contrasting measured versus intrinsic emission 
leads to a determination of the 7-7 opacity to VHE photons - hence, upon assuming a specific 
cosmology, we derive the EBL photon number density. Using, for each given source, different states 
of emission will only improve the accuracy of the proposed method. We demonstrate this method 
using recent simultaneous multi-frequency observations of the blazar PKS 2155-304 and discuss how 
similar observations can more accurately probe the EBL. 
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I. INTRODUCTION 

The Extragalactic Background Light (EBL), in 
both its level and degree of cosmic evolution, reflects 
the time integrated history of light production and 
re-processing in the Universe, hence the history of 
cosmological star-formation. Roughly speaking, its 
shape must reflect the two humps that characterize 
the spectral energy distributions (SEDs) of galaxies: 
one arising from starlight and peaking at A ~ 1 fjxa 
(optical background) , and one arising from warm dust 
emission and peaking at A ~ 100 /im (infrared back- 
ground) . 

Direct measurements of the EBL are hampered by 
the dominance of foreground emission (interplanetary 
dust and Galactic emission), hence the level of EBL 
emission is uncertain by a factor of several. 

One approach has been modeling the EBL aris- 
ing from an evolving population of galactic stellar 
populations: however, uncertainties in the assumed 
galaxy formation and evolution scenarios, stellar ini- 
tial mass function, and star formation rate have 
led to significant discrepancy among models (e.g., 
0, H HI |H, IH). These models have been used 
to correct observed VHE spectra and deduce (EBL 
model dependent) 'intrinsic' VHE 7-ray emissions. 

The opposite approach, of a more phenomenologi- 



cal kind, deduces upper limits on the level of EBL at- 
tenuation making basic assumptions on the intrinsic 
VHE 7-ray shape of AGN spectra: assuming, specif- 
ically, that the VHE photon index must be T > 1.5; 
e.g., (@, [H, El); but see ([H), or that the same- 
slope extrapolation of the observed Fermi /L AT HE 
spectrum into the VHE domain exceeds the intrin- 
sic VHE spectrum there ([§]). The only unquestion- 
able constraints on the EBL are model-independent 
lower limits based on galaxy counts (@, Q). It should 
be noted, however, that the EBL upper limits in the 
2-80/im obtained by [lj§] combining results from all 
known TeV blazar spectra (based on the assumption 
that the intrinsic T > 1.5) are only a factor «2-2.5 
above the absolute lower limits from source counts. 
So it would appear that there is little room for ad- 
ditional components like Pop III stars, unless we miss 
some fundamental aspects of blazar emission theory 
(which we never observed in local sources, however). 

An attempt to measure the EBL used the relatively 
faraway blazar 3C 279 as a background light source 
([24]), assumed that the intrinsic VHE spectrum was 
known from modeling and extrapolating the (histor- 
ical) average broad-band data. However, blazars are 
highly variable sources, so it's almost impossible to 
determine with confidence the intrinsic TeV spectrum 
- which itself can be variable. 

In this paper we propose a method to mea- 
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sure the EBL that improves on [24|] by making a 
more realistic assumption on the intrinsic TeV spec- 
trum. Simultaneous optical/X-ray/HE/VHE (i.e., 
eV/keV/GeV/TeV) data are crucial to this method, 
considering the strong and rapid variability displayed 
by most blazars. After reviewing features of EBL 
absorption (sect. 2) and of the SSC emission model 
(sect. 3), in sect. 4 we describe our technique, in sect. 5 
we apply it to recent multifrequency observations of 
PKS 2155-304 and determine the photon-photon opti- 
cal depth out to that source's redshift. In sect. 6 we 
discuss our results. 



II. EBL ABSORPTION 

The cross section for the reaction 77— >e ± is (12]), 
3 



c- 77 (-E,e) 



-a T (1-/3") x 



1G 

2/3(/3 2 -2) + (3-/3 4 ) In 



1-/3 



, (1) 



where ut is the Thompson cross section and 
fi=yj\ — (m e c 2 ) 2 /Ee. For demonstration purposes let 
us assume, following [24|, that n(e) is the local num- 
ber density of EBL photons having energy equal to 
e (no redshift evolution - as befits the relatively low 
redshifts accessible to IACTs), z e is the source red- 
shift, and Oo=l' the corresponding optical depth due 
to pair creation attenuation between the source and 
the Earth, is (see (HI) 
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III. BLAZAR SSC EMISSION 

In order to reduce the degrees of freedom, we use a 
simple one-zone SSC model (for details see [H, [29| ) . 
This has been shown to adequ ately describe broad- 
band SEDs of most blazars (|1C|12!|) an d, for a given 
blazars, both its ground and excited states ([27j]). 
The reason for the one-zone model to work is that in 
most blazars the temporal variability is clearly dom- 
inated by one characteristic timescale, which implies 
one dominant characteristic size of the emitting region 



The emission zone is supposed to be spherical with 
radius R, in motion with bulk Lorentz factor T at an 
angle 9 with respect to the line of sight. Special rela- 
tivistic effects are described by the relativistic Doppler 
factor, (5 = [T(l — /3cos#)] -1 . The energy spectrum 
of the emitting relativistic electrons is described by a 
smoothed broken power-law function of the electron 
Lorentz factor 7, with limits 71 and 72 and break at 
7 br . In calculating the SSC emission we use the full 
Klein-Nishina cross section. 

As detailed in [29], this simple model can be fully 
constrained by using simultaneous multifrequency ob- 
servations. Indeed, the total number of free param- 
eter of the model is reduced to 9: the 6 parame- 
ter specifying the electron energy distribution, plus 
the Doppler factor S, the size of the emission region 
R, and the magnetic field B. On the other hand, 
from observations ideally one can derive 9 observa- 
tional quantities: the slopes of the synchrotron bump 
after and above the peak a 1.2 (uniquely connected to 
^1,2)5 the synchrotron and SSC peak frequencies (24, c) 
and luminosities £ s ,Cj arL d the minimum variability 
timescale i V ar which provides an estimate of the size 
of the sources through R < ct var 5. 

Therefore, once the relevant observational quanti- 
ties are known, one can uniquely derive the set of SSC 
parameters. 



where x=(l—cos9), being the angle between the 
photons, and Ho is the Hubble constant. We further 
assume, again following [24|, that the local EBL spec- 
trum has a power-law form, n(e)oce -2 ' 55 . Then Eq.(l) 
yields t(E, z)otE 1 - 55 z^ with 77^1.5. 

This calculation, although it refers to an idealized 
and somewhat simplified situation, highlights an im- 
portant property of the VHE flux attenuation by the 
7vhe7ebl — ^e^e - interaction: r 77 depends both on 
the distance traveled by the VHE photon (hence on 
z) and on its (measured) energy E. So the spectrum 
measured at Earth is distorted with respect to the 
emitted spectrum. In detail, the expected VHE 7-ray 
flux at Earth will be: F(E)= (dl/dE) e - T -*-<^ (dif- 
ferential) and F{>E)=Jz (dl /dE') e- T ~<^ E ">dE' (in- 
tegral). 



IV. THE METHOD 

The method we are proposing stems from the con- 
sideration that both the EBL and the intrinsic VHE 7- 
ray spectra of background sources are fundamentally 
unknown. In order to measure the EBL at different 
z, one should single out a class of sources that is ho- 
mogeneous, i.e. it can be described by one same emis- 
sion model at all redshifts. This approach is meant to 
minimize biases that may possibly arise from system- 
atically different SED modelings adopted for different 
classes of sources at different distances. So we choose 
the class of source that both has the simplest emis- 
sion model and has the potentiality of being seen from 
large distances: blazars, i.e. the AGN whose relativis- 
tic jets point toward the observer so their luminosities 
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are boosted by a large factor and dominate the source 
flux with their SSC emission. Within blazars, we pro- 
pose to use the sub-class of "high-frequency peaked 
BL Lacs" (HBL), because their Compton peak can 
be more readily detected by IACTs than other types 
of blazar, and because their HE spectrum can be de- 
scribed as a single (unbroken) power law in photon 
energy unlikely other types of blazar ([ll.[l5l]). 

For a given blazar, our method relies on using, a 
simultaneous broad-band SED that samples the op- 
tical, X-ray, high-energy (HE: E > 100 MeV) 7-ray 
(from the Fermi telescope), and VHE 7-ray (from 
Cherenkov telescopes) bands. A given SED will be 
best-fitted, from optical through HE7-rays, with a 
Synchrotron Self-Compton (SSC) model. [Photons 
with E% 100 GeV are largely unaffected by EBL at- 
tenuation (for reasonable EBL models) as long as 
z^S 1.] Extrapolating such best-fitting SED model 
into the VHE regime, we shall assume it represents 
the blazar's intrinsic emission. Contrasting measured 
versus intrinsic emission yields a determination of 
e~ T ~t~t( E ' z \ the energy-dependent absorption of the 
VHE emission coming from a source located at red- 
shift z due to pair production with intervening EBL 
photons. Upon assumption of a specific cosmology, 
the final step is deriving the EBL photon number den- 
sity. 

Using, for each blazar, SEDs from different states of 
emission will improve the accuracy of the method by 
increasing the number of EBL measurements at that 
redshift. 



V. RESULTS: APPLICATION TO 
PKS 2155-304 

We apply the procedure described in Sect. 4 to 
the simultaneous SED data set of PKS 2155-304 de- 
scribed in 0- The data and resulting best-fit SSC 
model (from optical through HE7-rays) are shown 
in Fig. (JTJ) . The extrapolation of the model into the 
VHE 7-ray range clearly lies below the observational 
H.E.S.S. data, progressively so with increasing en- 
ergy. We attribute this effect to EBL attenuation, 
F hsi.E; z)—F cnl (E\ z) e~ T T~i( E ' z \ The corresponding 
values of t 77 (S; z) for £=0.23, 0.44, 0.88, 1.70 TeV 
and z—0.12 are, respectively, t 77 = 0.12, 0.48, 0.80, 
and 0.87 . 

We note that the SED analysis of Q was based on 
a slightly different SSC model, that involved a three- 
slope (as opposed to our two-slope) electron spec- 
trum. This difference may lead to a somewhat differ- 
ent decreasing wing of the modeled Compton hump, 
and hence to a systematic difference in the derived 
Tjj(E; z). That said, it's however interesting to note 
that the main parameters describing the plasma blob 
{B, 5, n c ) take on similar values in our best-fit analysis 
and in [2|. 

In Fig. ([2]) we compare our determination of 7Vv with 
some recent results ([8]) or upper limits ( 0,(3, HU). 
Whereas our values are generally compatible with pre- 
viously published constraints, we note that our val- 
ues closely agree with the corresponding values of [|| , 
which are derived from galaxy number counts and 
hence represent the light contributed by the stellar 
populations of galaxies prior to the epoch correspond- 
ing to source redshift z s - i.e., the minimum amount 
(i.e., the guaranteed level) of EBL. 



A. Best-fit procedure: \ 2 minimization 



In order to fit the observed optical, X-ray and HE 
7-ray flux with the SSC model, a % 2 minimization 
is used. We vary the SSC model's 9 parameters by 
small logarithmic steps. If the variability timescale 
of the flux, i var , is known, one can set R ~ ct v&I S, 
so the free parameters are reduced to 8. We assume 
here 7mi n =l: for a plasma with n o «O(10)cm~ 3 and 
£?«O(0.1) G (as generally appropriate for TeV blazar 
jets: e.g., [1, 0, El ) , this approximately corresponds 
to the energy below which Coulomb losses exceed the 
synchrotron losses (e.g., [13, Hl|) and hence the elec- 
tron spectrum bends over and no longer is power- 
law. However, in general 7 m ; n should be left to vary - 
e.g., cases of a "narrow" Compton component require 
7min>l (,30]). In order to reduce the run time of the 
code, the steps are adjusted in each run such that, a 
larger x 2 is followed by larger steps. 



VI. DISCUSSION 

The method for measuring the EBL we have pro- 
posed in this paper is admittedly model-dependent. 
However, its only requirement is that all the sources 
used as background beamlights should have one same 
emission model. In the application proposed here, we 
have used a one-zone SSC model where the electron 
spectrum was a (smoothed) double power law applied 
to the SED of the HBL blazar PKS 2155-304. While 
this choice was encouraged by the current observa- 
tional evidence fact that seem to HBLs have, with 
no exception, single-slope Fermi-L AT spectra ([HI]), 
we could have as well adopted the choice (01) of a 
triple power law electron spectrum in our search for 
the best-fit SSC model of PKS 2155-304's SED. Should 
the latter electron distribution be shown to provide a 
better fit to HBL Fermi-LAT spectra, then it would 
become our choice. In general, what matters to the 
application of this method, is that all source SEDs be 
fit with one same SSC model. 
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FIG. 1: Data (symbols: from [see:2]) and best-fit SSC model (solid curve) of the SED of PKS 2155-304. The best-fit 
SSC parameters are: n e = 150cm -3 , 7 br = 2.9 x 10 4 , 7 max = 8 x 10 s , cei = 1.8, a 2 = 3.8, R = 3.87 x 10 16 cm, 5 = 29.2, 
B — 0.056 G. The obtained values of R and S imply a variability timescale t V3lI ~ R/(c8), which is compatible with the 
observed value of ~12 hr. 



This work will be the subject of a forthcoming paper ( [17[ ) • 
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FIG. 2: Measured values of t 11 (E\ z) for _E=0.23, 0.44, 0.88, 1.70 TeV derived from comparing, for simultaneous 
observations of the HBL blazar PKS 2155-304 (z=0.12), the (EBL-affected) VHE7-ray data with the eV-through-GeV 
best-fitting SSC model extrapolated into the TeV domain. The curves represent, for redshifts z = 0.05, 0.1, 0.2, 0.3, 
0.4, and 0.5, the optical depth t 11 (E) according to [8| (top left), and upper limits to it according to [ll[ (top right), [TJ 
(bottom left), [l^ (bottom right). 



1520 (2008) 
[20] Y.Rephaeli, ApJ, 227, 364-369 (1979) 
[21] M.H. Salamon, & F.W. Stecker, ApJ, 493, 547-554 

(1998) 

[22] C.L. Sarazin, ApJ, 520, 529-547 (1999) 
[23] F.W. Stecker & O.deJager, A&A, 334, L85-L87 
(1998) 

[24] F.W. Stecker & O.C.deJager, & M.H. Salamon ApJ, 

390, L49-L52 (1992) 
[25] F.W. Stecker, M.A.Malkan & S.T.Scully, ApJ, 648, 

774-783 (2006) 



[26] F.W. Stecker, M.G. Baring & E.J. Summerlin, ApJ, 

667, L29-L32 (2007) 
[27] G. Tagliaferri, et al. (MAGIC collab.), ApJ, 679, 1029- 

1039 (2008) 

[28] F. Tavecchio & L. Maraschi, ApJ, 593, 667-675 (2003) 
[29] F. Tavecchio, L. Maraschi, & G. Ghisellini, ApJ, 509, 

608-619 (1998) 
[30] F. Tavecchio, G. Ghisellini, G. Ghirlanda, et al., MN- 

RAS, in press (2009) [arXiv: 0909.0651] 



eConf C091122 



